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Combined PET/computed tomography (CT) is of value in 
cancer diagnosis, follow-up, and treatment planning. For 
cancers located in the thorax or abdomen, the patient's 
breathing causes artifacts and errors in PET and CT 
images. Many different approaches for artifact avoidance or 
correction have been developed; most are based on gated 
acquisition and synchronization between the respiratory 
signal and PET acquisition. The respiratory signal is usually 
produced by an external sensor that tracks a physiological 
characteristic related to the patient's breathing. Respiratory 
gating is a compensation technique in which time or 
amplitude binning is used to exclude the motion in 
reconstructed PET images. Although this technique is 
performed in routine clinical practice, it fails to adequately 
correct for respiratory motion because each gate can mix 
several tissue positions. Researchers have suggested 
either selecting PET events from gated acquisitions or 
performing several PET acquisitions (corresponding to a 
breath-hold CT position). However, the PET acquisition time 
must be increased if adequate counting statistics are to be 
obtained in the different gates after binning. Hence, other 



researchers have assessed correction techniques that take 
account of all the counting statistics (without increasing 
the acquisition duration) and integrate motion information 
before, during, or after the reconstruction process. Here, 
we provide an overview of how motion is managed to 
overcome respiratory motion in PET/CT images. Nucl Med 
Commun 35:113-122 © 2014 Wolters Kluwer Health | 
Lippincott Williams & Wilkins. 
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Introduction 

PET imaging has proven value in cancer diagnosis, follow- 
up, and patient management [1,2]. It can also be used to 
design and adjust treatment plans in radiotherapy and 
chemotherapy [3]. 

The introduction of hybrid PET/computed tomography 
(CT) scanners over the last decade has improved the 
accuracy of attenuation correction and the ability to 
localize organs and tumors [4,5]. However, a number of 
physical and physiological phenomena can affect the 
accuracy of PET and CT images. Indeed, it has been 
shown that the diaphragm can move by as much as 20 mm 
in the craniocaudal axis [6] and that the liver can move by 
an average of 11mm [7] during quiet breathing in the 
supine position. It is also well known that breathing 
induces rotational and/or translation movements in 
thoracic and abdominal organs to a varying extent [8-10]. 
For example, upper areas of the lungs are less subject to 
motion compared with the lower parts [6]. Moreover, 
Rodarte et al. [8] showed that even though the middle 
lobe of the canine lung and the lingula extend toward the 



diaphragm these structures are less subject to respiratory 
motion compared with the lower lobes. 

CT images are acquired within a few seconds, whereas 
PET images are acquired over several minutes for each 
axial field of view. Even though the various CT slices in 
a free-breathing CT volume can be considered to be 
nonsynchronized 'snapshots', they will contain tissues 
from different respiratory states because the timeframe of 
a CT acquisition is roughly the same as that of the 
respiratory cycle (RC). Inversely, PET images of a moving 
uptake represent motion averaged over a few minutes. 
This difference in timeframe may induce regis- 
tration errors, respiratory motion CT artifacts, and/or 
blurred PET images [11]. Furthermore, PET smearing 
effects and the attenuation correction derived from CT 
may in turn induce erroneous quantification of the 
standardized uptake value in PET. 

Several procedures and devices have been developed with a 
view to avoiding these problems. Here, we review current 
methods for dealing with respiratory motion in PET/CT 
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Respiratory motion tracking 

Tracking techniques rely on an external device to estimate 
respiratory motion during PET/CT examinations. There is 
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a proven correlation between respiratory motion and the 
displacement of internal organs [12-14]. Respiratory 
motion can be estimated by monitoring a physiological 
characteristic (such as the displacement of the thoracic 
cage or the volume of air exhaled) or a physical parameter 
(such as temperature or displacement of the skin surface). 
For analog sensors, the signal must be converted into a 
digital signal, which will be referred to as the respiratory 
signal (RS) throughout this review. Gated PET acquisition 
involves synchronization of the PET acquisition with the 
RS. 

Signal measurement 

Several solutions have been implemented on commercial 
imaging systems by the PET/GT device manufacturers, 
whereas academic researchers have assessed a number of 
other techniques. 

The AZ-733V (Anzai Medical Corp., Tokyo, Japan) is a 
pressure sensor integrated into an elastic chest belt 
(Fig. la). The pressure sensor can thus detect the 
displacement of the abdominal wall as a function of 
pressure variations during the RC (i.e. low pressure 
during expiration and high pressure during inspiration). 
An analog-to-digital converter yields a signal that is 
recorded as a text file and can be used for offline 
processing [16]. 

Another widely used method involves the Real-Time 
Position Management (RPM) system (Varian Medical 
Systems, Palo Alto, California, USA). Two infrared 
reflective markers are placed on a plastic box positioned 
on the patient's thorax. A video camera placed at the end 
of the bed tracks the markers' motion during the RG 
(Fig. lb). The markers' displacement is recorded by the 
video camera and yields an RS for further processing [16]. 

Otani era/. [18] compared the AZ-733Vand RPM external 
respiratory tracking devices and found that the respective 
RSs were well correlated. 

Sensors measuring air temperature changes during 
respiration can also be used. Indeed, the air temperature 
in the upper airways during inspiration is lower than that 
during expiration because the air is warmed during its 
passage through the lungs. A high time resolution is 
needed for this type of sensor. In one embodiment of this 
system, a high-sensitivity thermistor can be mounted 
inside a conventional oxygen mask. In another variant, a 
probe is placed close to the patient's nostrils (BioVet GT1 
System; Spin Systems, Brisbane, Australia) (Fig. lc) [16]. 
Boucher et al. [19] have demonstrated the clinical 
feasibility of this type of device. However, this tracking 
device has not yet been implemented on commercial 
PET scanners and has not been greatly characterized in 
the literature. 

A spirometer can be used to estimate the volume of air 
inhaled or exhaled during breathing. It can be placed 



close to the patient's nostrils (PMM Spirometer; Siemens 
Medical Systems, Erlangen, Germany [16]) or mouth 
(GPX Spirometer; Medgraphics, St Paul, Minnesota, 
USA [17]) (Fig. Id). However, use of this system with 
a clinical PET/GT gantry has not been investigated yet. 
The main drawbacks of this system are discomfort for the 
patient and the inability to provide simultaneous oxygen 
assistance. Similarly, Didierlaurent et al. [20] connected a 
pneumotachograph (which measures the change in lung 
volume over time via the exhaled respiratory flow) to a 
PET/CT gantry. 

Synchronizing a respiratory signal with PET or 
computed tomography data 

The RS can be synchronized with four-dimensional (4D) 
PET or CT data. In list-mode (LM) acquisitions, 
coincidence events and time tags are recorded in real 
time. As shown by Bruyant et al. [21], a tag can be stored 
within the LM file whenever the user-defined RS 
threshold is reached. Given that the LM retains time 
data, PET synchronization and processing can be 
performed after the acquisition. 

Compensation techniques 

The respiratory gating methods based on motion tracking 
mostly use LM PET data to take advantage of the 
temporal component. There are two ways to process 
respiratory-gated PET data: dividing the RS into multiple 
phases (in multiphase or multibin methods) or consider- 
ing only a single phase (in single-bin methods). 

Multibin methods 

The RS can be processed as a function of time or 
amplitude. In both time-based and amplitude-based 
methods, the RG can be divided into equal-sized or 
different-sized bins. We shall detail each of these 
approaches in the following sections. 

Time-based methods 

In the time-based approach, each RG (as delimited by 
two gating tags in the LM) is divided into several 
bins [22-27]. Several different approaches have been 
suggested. Dawood et al. [28] set a fixed time for all bins 
and for all RCs (Fig. 2a). However, depending on the RG 
length, some PET data may be lost. Manufacturers have 
offered time-based, respiratory-gated PET processing 
based on the division of each RG into the same number 
of bins (Fig. 2b) [29]. In this method, users have to 
define a range of acceptable RG frequencies. Data from 
abnormally short or long RGs are excluded from the 
analysis. 

The number of bins is an important determinant of image 
quality. In fact, there is a tradeoff between residual 
motion and counting statistics: the lower the number of 
bins, the weaker the blurring removal but the higher 
the signal-to-noise ratio in each bin. Conversely, a high 
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Fig. 1 




(a) The AZ-733V pressure belt (Anzai Medical Corp.). (b) The Real-Time Position Management (RPM) system (Varian Medical Systems) (reprinted 
with permission from AAPM and Dr Sadek Nehmeh [1 5]). Copyright AAPM, College Park, MD, USA. All permission requests for this image should be 
made to the copyright holder, (c) The BioVet CT1 System (Spin Systems) (reprinted with permission from Dr Axel Martinez-Moller [16]). Copyright 
Technical University Munich, Germany. All permission requests for this image should be made to the copyright holder, (d) The CPX Spirometer 
(Medgraphics) (reprinted with permission from AAPM and Dr Bhudatt Paliwal [17]). Copyright AAPM, College Park, MD, USA. All permission 
requests for this image should be made to the copyright holder. 



number of bins translates into negligible residual motion 
and a marked decrease in the signal-to-noise ratio. On the 
basis of Monte Carlo methods, Vauclin et al. [30] 
suggested that, although the accuracy in lesion volume 
estimation increases with the number of bins considered, 
there is a corresponding decrease in overall lesion 
detectability because the total count statistic was divided 
by the number of gated bins. Dawood et al [31] reported 
that six respiratory bins constituted the optimal tradeoff 
between residual motion in an individual bin on one hand 
and image quality on the other. 

Furthermore, the degree of motion compensation differs 
from one bin to another. This is due to differences in 
tissue velocity in the different phases of the RG; some 
bins correspond to a very steep part of the RS, where the 
tissue velocity is high. Next, uptake restitution of the 
moving target is flawed by high residual motion. More- 
over, respiratory amplitude varies significantly from one 
individual to another. Hence, events from different tissue 
positions may become mixed within a single bin, leading 
to a blurred image of the target. To minimize these 
differences and regularize the RG, some systems offer 
'vocal coaching'. This method is only partially effective, 
as 'coached' respiration is not always regular [15,32]. 

Many researchers have compared gated PET images with 
ungated PET images (i.e. those acquired in the absence 
of RG synchronization) [22-27]. Table 1 summarizes the 



results associated with time-based methods [24-27]. 
Depending on the lesion size, respiratory motion can lead 
to major errors in uptake detection. For motionless 
lesions, the detection limit is roughly two or three times 
the tomograph's spatial resolution [33]. The quantifica- 
tion bias increases as the lesion's largest axis de- 
creases [34]. Adding a motion component raises the 
limit of detection and increases quantification bias. 

The time-based approach has been implemented by all 
manufacturers because it does not require access to the 
RS. Indeed, the RS is generally obtained by third-party 
devices and is neither available nor processed on 
reconstruction workstations. 



Amplitude-based methods 

In this approach, the RS is divided into several bins as 
a function of its magnitude. Minimum and maximum 
thresholds define the range of motion (i.e. the range of 
amplitudes) to be processed. In one method, each bin 
contains equal ranges of amplitude [28,29]. Use of an 
amplitude-based method means that the various bins are 
unlikely to contain the same number of coincidences 
(Fig. 2c). Indeed, the number of statistics in a steep part 
of the RS will be lower than that in the end-expiratory 
phase. To overcome this issue, variable bin processing has 
been suggested; each bin contains the same statistics but 
with variable remaining motion amplitude [28]. 
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Fig. 2 
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Table 1 Overview of studies using time-based methods 





Werner et al. [24] 


Lupi et al. [25] 


Garcia Vicente et al. [26] 


Kasuya et al. [27] 


Gating 












Tracking device 


RPM 


RPM 




RPM 


AZ-733V 


Number of bins 


8 


6 




6 


5 


Lesions 












Site 


Lung 


UL LL 


UL 


LL 


Pancreatic head Pancreatic groove 


N 


18 


19 3 


29 


13 


13 2 


suv max 












Gated 


11.8 (±5.5) 


15.7 (±13.3) 10.1 (±10.9) 


2.2 (±0.8) 


2.5 (±0.9) 


5.8 (±4.3) 


Ungated 


9.2 (±4.8) 


9.6 (±7.1) 5.9 (±6.0) 


1.3 (±0.6) 


1.4 (±0.4) 


5.4 (±3.9) 



LL, lower lung; RPM, Real-Time Position Management system; SUV max , maximum standardized uptake value; UL, upper lung. 



As for the time-based methods, the number of bins is 
an important determinant of image quality. Dawood 
et al. [31] showed that the optimal configuration 
contained eight respiratory bins. However, Bettinardi 
et al. [35] emphasized that the number of gates needed is 
related to a lesion's size and displacement. 

Although the amplitude-based approach is superior to the 
time-based approach [28,29], the workstations require 
access to the RS during data processing. Multibin 
approaches can describe a lesion's displacement during 
the RC, which can be of value for radiotherapy 
planning [36]. However, these two approaches do not 
provide adequate attenuation correction. Indeed, GT 
images, usually acquired in free-breathing mode in a 
single volume, show a number of respiratory artifacts and, 
furthermore, do not correspond to a defined respiratory 
state. These two phenomena introduce errors into the 
attenuation correction and quantification of the up- 
take [29]. 

Single-bin methods 

The deep-inspiration breath-hold method 

Nehmeh etal. [37] and Meirelles et al. [38] developed the 
deep-inspiration breath-hold (DIBH) method. A device is 
used to track RS and deliver voice instructions (e.g. 
'breathe in', 'hold', and 'relax'). A 15-s GT scan chest is 
acquired while the patient holds his/her breath in deep 
end-inspiration. The PET stage consists of nine inde- 
pendent, 20-s breath-hold acquisitions (again during 
deep end-inspiration). The patient is allowed to relax 
for 20 s between each acquisition. Finally, these acquisi- 
tions are summed, corrected for attenuation (using the 
GT data), and reconstructed with a conventional algo- 
rithm (Fig. 3). 

With the DIBH method, both PET and CT data are 
acquired at the same tissue position during a deep end- 
inspiration breath-hold; this guarantees accurate attenua- 
tion correction. However, repeated 20-s acquisitions may 
be stressful and tiring in dyspneic patients. Even when 
breathing is coached and monitored, there is no guarantee 
that the patient will hold his/her breath at the same point 
in all the acquisitions. To avoid the need for repeated 
breath-hold acquisition, two groups have shown that 
it is feasible to acquire PET data during a single 



DIBH [39-41]. Kawano et al. [39] and Torizuka 
et al. [41] respectively considered that a breath-hold of 
30 or 20 s is sufficient. However, the body phantom 
studies used in these two studies do not fully model 
clinical reality because of high concentration of activity in 
the spheres and the absence of noise in the thoracic 
cavity. 

The DIBH method has been assessed in routine clinical 
practice (mainly for lung imaging) [38-43]. Indeed, the 
method allows the fully expanded lungs to be imaged and 
small tumors can be detected easily. Other researchers 
have evaluated this method for imaging the abdo- 
men [43,44]. The main findings of these DIBH studies 
are summarized in Table 2. 

The computed tomography-based method 

Daouk etal. [45] and Fin et al. [46] have described a GT- 
based method, in which a 10-min LM PET acquisition is 
performed. The patients then hold their breath at the 
normal end-expiration point for 6 s, during which a 
breath-hold GT (BH-GT) scan is performed. The RS is 
recorded throughout the two sessions. The researchers 
chose this normal end-expiration because it is the most 
reproducible phase of the RC [47]. Given that end- 
expiration lasts longer than any other phase, selecting the 
events recorded during this phase improves counting 
statistics and provides the best image quality. In end- 
expiration, a plateau on the RS corresponds to a specific 
tissue position. An event selection range (ESR) is defined 
around this plateau. The portions of the RS falling within 
this ESR thus define the portions of the LM that should 
be conserved. The selected PET data are attenuation- 
corrected with the BH-GT data during reconstruction. 
Placing the ESR around the breath-hold position defined 
by BH-GT acquisition provides accurate matching 
between anatomical and functional images and yields 
accurate attenuation correction (Fig. 4). 

The drawback of using end-expiration is the low lung 
volume in this phase of the RG, which may make it more 
difficult to detect small tumors. Furthermore, the GT- 
based method rejects a large proportion of the acquired 
PET data. Hence, a 10-min acquisition is needed to 
achieve the same quality as a standard 3-min acquisi- 
tion [48]. 
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Table 2 Overview of studies using the deep-inspiration breath-hold method 



Meirelles 
et al. [38] 



Kawano 
et al. [39] 



Kawano 
et al. [40] 



Torizuka et al. [41] 



Daisaki 
et al. [42] 



Nagamachi et al. [43] 



Nagamachi 
et al. [44] 



Protocol 

Number of 9 
breath-holds 

Time (s) 1 80 

Lesions 

Site Lung 

N 62 
Diameter (mm) 2-96 a 



suv max 

Gated 
Ungated 



8.0 (NA) 
7.3 (NA) 



1 

30-125 
UL LL 



1 

30-143 
UL LL 



78 30 40 55 

12-40 a 11-40 a 16-60 a 6-72 a 



10.5 7.8 9.4 10.6 

(±6.2) (±3.9) (±5.3) (±4.9) 

9.1 5.3 7.5 6.8 

(±5.5) (±2.7) (±4.7) (±3.6) 



1 

20 



UL 



LL 



22 25 

28.8 36.3 
(±14.8) (±15.5) 



8.6 
(±4.8) 

7.0 
(±3.9) 



13.6 
(±8.3) 

9.3 
(±4.9) 



>2 

120 
Lung 
32 

25.9 (±6.9) 

4.1 (±3.1) 
4.3 (±2.8) 



Esophagus, 
chest, lung 

NA 

11.7 (±8.5) 



4 

120 

Liver, bile duct, 

pancreas 
NA 



12.0 (±8.9) 
9.1 (±7.5) 



13.5 (±7.3) 
9.5 (±4.4) 



4 
120 

Liver Pancreas 



NA 



18.4 

(±4.2) 
11.6 
(±3.0) 



NA 



6.9 
(±1.5) 

6.0 
(±1.3) 



LL, lower lung; NA, not available; SUV max , maximum standardized uptake value; UL, upper lung. 
a Data provided with the range of results. 



When applied to the lungs [46,49] and the liver [50,51], 
the GT-based method was found to be more sensitive 
than ungated imaging. Table 3 summarizes the differ- 
ences in maximum standardized uptake values in these 
studies. 

Chang et al [52] have proposed a solution in which 
respiratory motion amplitude is captured during a whole- 
body GT acquisition. The RS amplitude is picked up at 
the entry and exit of the bed step where the lesion is 
localized in order to define an amplitude range (AR). 
Unlike the previously mentioned method, a gated, 
triggered LM acquisition then gathers only PET events 
corresponding to the defined AR. These PET data are 
used to generate a single motion-free volume. The 
drawback of this solution relates to its dependence on 
the patient's breathing pattern; if the patient breathes 
quickly, the influence of respiratory motion cannot be 



fully excluded. Moreover, the patient may not breathe in 
the same way for the GT and PET acquisitions. In such 
an event, the AR determined in the GT scan will not 
correspond to the respiratory amplitudes at the lesion's 
site according to the PET acquisition. 

Data-driven approach 

Nehmeh et al [53] developed a gating method using a 
single tissue position. The tissue position is determined 
with respect to an external 18 F-FDG point source that 
appears on each image of a dynamic PET acquisition. On 
the first reconstructed frame, the point source's position 
is defined as the reference, and a region of interest (ROI) 
is drawn around it. The ROI is then carried over onto all 
frames but only those in which the point source falls 
within the ROI are considered. Finally, the corresponding 
PET raw data are summed and reconstructed. 
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Table 3 Overview of studies using the computed tomography-based method 





Fin et al. [46] 


Daouk et al. [49] 


Daouk et al. [50] 


Fin etal. [51] 


Lesions 
















Site 


UL LL 




Lung 




Liver 




Liver 


N 


5 8 




48 




74 


4 


27 


Diameter (mm) 


14.0 (±5.3) 


< 15 




>15 


15-25 


10-14 


> 15 


suv max 
















Gated 


4.9 (±2.8) 


2.3 (±1.8) 




5.7 (±2.5) 


8.5 (±3.5) 


5.3 (±0.9) 


7.7 (±3.2) 


Ungated 


3.8 (±2.3) 


1.9 (±1.3) 




5.2 (±2.5) 


7.6 (±3.1) 


4.8 (±0.5) 


7.1 (±2.5) 



LL, lower lung; SUV max , maximum standardized uptake value; UL, upper lung. 

Correction techniques 

To avoid increasing the duration of PET acquisitions (due 
to the loss of statistics related to compensation techni- 
ques), some researchers have developed techniques to 
correct for respiratory motion. To the best of our 
knowledge, these techniques have not yet been used in 
routine clinical practice. Motion correction is usually 
performed by introducing a deformation field before, 
during, or after the tomographic reconstruction. Deter- 
mination of the deformation field is beyond the scope of 
this paper but we shall detail these three approaches 
(before, during, or after reconstruction) in the following 
sections. 

Prereconstruction correction methods: data-driven 
approaches 

Prereconstruction correction involves changing the posi- 
tion of lines of response (LORs) [54-57] as a function of 
a precomputed deformation field. However, the LORs are 
displaced by means of affine transformations only; this 
means that the technique is well suited to the correction 
of head motion [58,59] but not to the correction of the 
nonrigid motion typically seen during breathing. The loss 
of some LORs may occur during repositioning (e.g. 
outside the gantry or between the detectors in 3D 
acquisitions) [57]. Moreover, the exact projection of the 



deformation field in raw space data may require a 
rounding step in LOR repositioning, which introduced 
noise into the final raw data [56]. 

Motion correction during reconstruction 

With the availability of fast computers, PET images are 
usually obtained by applying expectation maximization 
algorithms [60-63]. Continuous-object reconstructions 
are generally based on discrete-data formulations through 
parallelepiped functions (i.e. voxels) and a system matrix, 
which links acquired data and the object's basis functions. 

To correct for respiration motion, several groups have 
suggested the computation of a system matrix for each 
respiratory state [64-67]. The main drawback of these 
approaches relates to the rigid nature of the basis 
function commonly used. Indeed, voxels do not deform 
along with the motion correction, which generates 
discontinuities between voxels. The solution presented 
by Reyes and colleagues [68,69] uses spherical basis 
functions inscribed in voxels. Basis functions could be 
approximated by deformable ellipsoids, which avoid the 
discontinuities caused by tissue deformation. 

Another way of correcting for respiratory motion is to take 
account of a-priori knowledge about an object (edges, 
anatomical features, etc.) by incorporating spatiotemporal 
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regularization into the reconstruction algorithm [70-74]. 
Grotus et al. [75] integrated temporal regularization with 
a temporally adaptive basis function. However, one 
drawback of these solutions relates to the determination 
of scalar weighting parameters, which can be performed 
empirically on phantoms but is not necessarily suited to 
clinical situations. 

Postreconstruction motion correction 

This method relies on summing gated images (obtained 
with a multibin approach) after a registration step 
[76-81]. In principle, this is an attractive approach, 
although PET image registration is complicated by the 
spatial resolution and by partial volume effects. In 
addition, this method suffers from the limitations 
described above in the section on multibin compensation. 
Moreover, the summed image is affected by the sum of 
the noise levels in each image. Nevertheless, registration 
of the gated PET images enables the determination of 
the deformation field required by the above-mentioned 
motion correction methods. 

Future trends and conclusion 

At the outset, respiratory motion compensation and 
correction techniques were aimed at improving diagnosis 
and follow-up during nuclear medicine examinations of 
the thorax and abdomen. However, recent advances in 
radiotherapy mean that metabolic images are taken into 
account in treatment planning [82]. Moreover, intensity- 
modulated radiotherapy planning and dose painting 
require information on the respiratory motion of lesions. 
In this context, multibin compensation techniques are 
the only methods capable of describing the lesion's path 
during the RC. Nevertheless, specific attenuation maps 
for each phase are required to describe intrinsic hetero- 
geneities within lesions. One solution is provided by 4D 
GT, albeit at the cost of a substantial increase in the 
radiation dose and acquisition time. Now that combined 
PET/MRI scanners are becoming more widely available, a 
few reports have emphasized the value of MRI for correct 
PET images for respiratory motion [83-85]. In MRI, the 
patient is not exposed to ionizing radiation. Moreover, the 
fact that most recent gantries perform PET and MRI 
acquisitions simultaneously (rather than sequentially, as 
in PET/GT scanners) improves the quality of registration. 
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